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Abstract 



& ' We consider a millicharged neutrino with nonzero magnetic moment in 

^ the presence of rotating and magnetized background matter. The exact so- 

^-i lution of the corresponding modified Dirac equation for the neutrino wave 

function is found. The neutrino energy spectrum is obtained and the effect 
of neutrino energy quantization is discussed in details. We introduce a new 
kind of spin operator which is a superposition of longitudinal and transverse 
polarizations operators for description of the neutrino spin properties in the 
considered background environment. Within the quasi-classical approach to 
the problem, radius of the neutrino orbits is derived and the effective "matter 
induced Lorentz force" is introduced. It is shown that in the considered envi- 
ronment, and also in matter with nonzero gradient of density, neutrino moves 
with acceleration. In this case a new type of the electromagnetic neutrino 
radiation (termed "light of millicharged neutrino", LCv) can be produced. 
The considered problem is of interest for astrophysical applications, in par- 
ticular in connection with the recently reported hints of ultra-high energy 
neutrinos E — 1 10 PeV observed by IceCube jlj]. 
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1. Introduction 



Neutrino electromagnetic properties are among the most intriguing and 
exciting issues of modern particle physics. Within the Standard Model neu- 
trinos are massless and have "zero" electromagnetic properties. However, as 
it is well known in different extensions of the Standard Model a massive neu- 
trino has non trivial electromagnetic properties (for a review of the neutrino 
electromagnetic properties see 0, l^]). That is why it is often claimed that 
neutrino electromagnetic properties open "a window to new physics" [i| • 

A neutrino magnetic moment, as expected in the easiest generalization 
of the Standard Model, is very small and proportional to the neutrino mass, 
fi v ~ 3 x 10~ 19 fi B (m u /l eV) with fi B being the electron Bohr magneton. 
Much greater values are predicted in various other Standard Model general- 
izations (for details see [3,0]). 

It is usually believed that the neutrino has zero electric charge. This 
can be attributed to gauge invariance and anomaly cancellation constraints 
imposed in the Standard Model. However, if the neutrino has a mass, the 
statement that a neutrino electric charge is zero is not so evident as it meets 
the eye. In theoretical models with the absence of hypercharge quantization 
the electric charge also gets "dequantized" 0-0] and as a result neutrinos 
may become electrically millicharged particles. Charge dequantization for 
a massive Dirac neutrino may occur, for instance, in the Standard Model 
extensions with right-handed neutrinos vr and also in a wide class of models 
that contain an explicit U(l) symmetry. It should be mentioned that only 
Dirac neutrino may have electric charge. 

The most severe experimental constraints on the electric charge 0, 0] and 
the magnetic moment , [ll| of the neutrino are given by 



where q u is a module of the neutrino millicharge and fis = eo/(2m e ). 

Due to obvious "smallness" of neutrino electromagnetic properties it is 
reasonable to expect that these properties manifest themselves most spec- 
tacular under the influence of extreme external conditions, namely in strong 
electromagnetic fields and dense background matter. These conditions can 
be easily found in astrophysics and cosmology settings. 

In this letter we consider new possible phenomena of neutrino electro- 
magnetic properties manifestation, in particular effects of nonzero magnetic 




q v < l(r 21 e , fi„ < 2.9 x 10~ n fi B , 
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moment and electric milllicharge, in case of the particle motion in dense 
magnetized matter. 

It is well known 12Hl5| that in case of particles motion in external high 
intensive electromagnetic fields the most adequate way to account for fields 
influence is based on the method of exact solutions of quantum equations for 
particle wave function. The same method of exact solutions can be used (see, 
for instance, [16} ) to describe particle (neutrinos and electrons, in particular) 
motion in dense background matter. In connection with evaluation of quan- 
tum theory 17, li[ of spin light of neutrino in matter SLu [l9| the method of 
exact solutions was applied for describing neutrino quantum states in dense 
matter. 

Most pronouncedly the power of the method of exact solutions is demon- 
strated in studies of different particle interactions in extreme conditions pecu- 



liar for astrophysical and cosmological settings [20H24| . Most recently exact 



solutions for a neutrino (and for an electron) in presence of matter and elec- 
tromagnetic fields of different configurations had been evaluated in 25M27 



see also [lJrM |28|) 



In this letter we consider for the first time a new form of the modified 
Dirac equation that describes the motion of an electrically millicharged neu- 
trino with anomalous magnetic moment (or any charged fermion including 
an electron) under the influence of an external electromagnetic field and also 
accounting for weak interaction with the background matter. In the next 
section we discuss in details derivation of the corresponding exact solution 
for the particle wave function. The exact expression for the energy spectrum 
is also obtained. Note that in order to treat the particles spin properties we 
introduce a new type of spin operator. 

In Section 3 we consider more specific background matter conditions and 
derive the exact solution for the millichaged neutrino wave function in case 
of magnetized and rotating matter. These particular external conditions are 
typical for the astrophysical environment that can be found, for instance, 
in neutron stars. An interesting new result is the neutrino energy spectrum 
quantization, the effect originated due to both electromagnetic and weak 
interactions of the neutrino with the background. 

In Section 4 within the quasi-classical interpretation using the neutrino 
energy spectrum we introduce the "effective Lorentz force" that defines the 
neutrino motion. This new force is induced not only by the electromagnetic 
interaction but by the weak interactions with rotating matter as well. Matter 
induced "charge", "electric" and "magnetic" fields are introduced. 
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In Section 5 a new possible mechanism of the electromagnetic neutrino 
radiation in nonuniform matter due to neutrino millicharge ( "Light of (milli) 
Charged Neutrino", LCv) is discussed. In particular, it is shown that the 
LCv can be emitted in the case of ultra-high energy millicharged neutrino 
moving in extreme astrophysical media with nonzero gradient of matter den- 
sity. This phenomenon can occur in extremely high density matter mostly 
composed of neutrons, such as neutron stars. The influence of the electron 
matter component on the LCv power of different neutrino flavors is also 
discussed. 

2. Millicharged neutrino wave function in matter and magnetic 



Consider a millicharged neutrino with nonzero magnetic moment moving 
in dense matter composed of neutrons in the presence of external magnetic 
field. In the most general case the modified Dirac equation for a fermion 
wave function with nonzero anomalous magnetic moment exactly accounting 
for interaction with matter and electromagnetic field is 



where F^ v is the electromagnetic field tensor and /j, is the particle's anomalous 
magnetic moment. For definiteness we suppose that a fermion has a negative 
charge q = —q . 

The matter effective potential V m = |7 m (q + 75)/^ includes both the 
neutral and charged current interactions of the fermion with the background 
particles. The explicit form of / M depends on matter composition and the 
background particles densities, speeds and polarizations. For the case of 
unpolarized matter at rest one obtains 



where Gp is the Fermi constant. In equations (j2j) and flSj) one should use 
charge, mass and magnetic moment of the considered fermion. For the 
electron neutrino q = c Ve = 1 and n = —n n , whereas for the electron 
q = c e = 1 — 4 sin 2 9w and n = n n . Here n n is neutron density. 



field 



lM + qoA") - \l,{ Cl + l5 )F ~ \^^ v - m ) = 0, (2) 




(3) 
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Equation (j2J) can be written in the Hamiltonian form i-^>{x) = H^/(x) 
with the Hamilton operator given by 



H = 7o7(p + QoM x )) + 7o™ + lo^ixB 



2^2 



ci + 75) 



(4) 



where p = i(d x , d y , d z ) and the magnetic filed potential is = (0, — 0). 

Using the standard representation of the 7-matrixes the Hamiltonian (jl]) 
can be re-written in the matrix form 



H 



(Gn + m + iiB A - ^ 

Gn + m- fiB A + iP 2 - P, - ^ 

P. - % P\-iP% Gn-m-fiB 

-P 3 - Gn-m + fiBj 



3 2 

V A + *A 



A - *A \ 

Gn 





(5) 

where P = p + qoA(x) is the "extended" momentum operator, and we use 
the notations 

G = £|, G=|g. (6) 

In the coordinate representation (the cylindrical coordinate system is used) 
the operators are given by the following expressions: A ±«A — ~i> e±%ip {jh =t 
"F lf r )> A — — The eigenvectors of the Hamiltonian (jl]), which 
satisfy the equation H^(x) = po^/(x), can be written in the form 



V(x) = Ce- i(Pot - p3z) 



tC 2 C l s (zfr 2 )e il * 

C 3 £i- 1 (T r2 ) e<(i " 1)v 
V iC 4 £ l s (3f-r 2 )e il v ) 



(7) 



where £i(%-r 



0,1,2...) and C is the 



, , , is a Laguerre function (JV = / + s 
normalization constant. The indeterminacy in the solution ([7]) due to the 
presence of undetermined arbitrary coefficients (i = 1, 2, 3, 4) still remains 
because the particle spin properties up to this point have not been specified. 
Using the Hamiltonian fl5]) and the solution of eq. ([2]) in the form given by 
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(J7j) we can get the following system of equations for the coefficients 



((Gn - po) + (m + d + (p 3 - C 3 + y/2Nq BC 4 = 0, 

((Gn - po) + (m - /xfl)) G 2 - fp 3 + ^ J G 4 + y/2Nq BC 3 = 0, 

((Gn - po) - (m+ fiB)j G 3 + I p 3 - j G x + y f 2N^BC 2 = 0, 

((Gn - po) - (m - G 4 - (^p 3 + G 2 + v/^iV^BG! = 0. 

To advance in evaluation of the exact solution of eq. (|2]) one should de- 
scribe the particle spin properties and to construct the spin operator that 
commutes with the Hamiltonian. For the considered case of both the mag- 
netic field and background matter presence this problem has never been 
studied before. To solve this problem we recall the known spin operators 
for the case of the particle motion in a pure magnetic field (without matter 
presence) and in matter (without a magnetic field) and also account for the 
fact that the new spin operator we are searching for should turn into one of 
the known spin operators in the corresponding limits of the magnetic field 
or background matter is "swiching off" : fiB — > or Gn — > 0. Thus we arrive 
to the following new spin operator 

Gn 

S = sin a S tT — cos aSi ong , sin a = — -. 2 , (9) 

which is a weighted superposition of the operators of the longitudinal and 
transverse polarizations, 

S tr = — , S long = Z 3 + -(° 7 )[<7XP] 3 , (10) 

m m \ctq u / 

where are the Pauli matrixes. The value of angle a is fixed by the values 
of the effective matter density n and the strength of the magnetic field B. 

It is easy to check that the introduced new spin operator (Q commutates 
with the Hamiltonian (j3J). Thus, these two operators have the same set of 
the eigenvectors: S^(x) = S^(x). 
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We obtain the system for the coefficients Cj which describes spin proper- 
ties of the wave function (171): 



(m cos a — P3 sin a + mS)Ci + \f2NqoB(C4, cos a — C% sin a) = 0, 
(m cos a — P3 sin a + mS)C 3 — \/2Nq B(C2 cos a + C\ sin a) = 0, 
[m cos a — p 3 sin a — mS)C 2 + \/ 2Nq B (C3 cos a + Ci sin a) =0, 
(m cos a — p 3 sin a — mS)C4 — \f2Nq$B(Ci cos a — C 3 sin a) = 0. 



Demanding for non-trivial solution of this system we obtain the eigenvalues 
of the spin operator ©, 



S = C — \/ (m cos a — P3 sin a) 2 + 2Nq B, C = ±1. 



(12) 



After that we can solve the systems of equations (|S]) and (TTTT) and get the 
energy spectrum 



Po = z 



i 



V\ + 2Nq B + m 2 + ( ^ + (/i-B) 2 - 2mSy + (//5) 2 +Gn 

(13) 



and the exact expressions for the spin coefficients 



1 / mcosa -p 3 sina , 

Cl = „ V 1 c V 1 + S1I1 (« + P)> 

2 V mo 



1 . , L m cos a — p 3 sin a n —, 7^ 

C 2 = -5iC\/ 1 + ^ Vl + sin , 

2 V mi 



C3 = -&>\/ 1 



m cos a — p 3 sin a 
mS 



'I - sin(a + p), 



(14) 



1.^1 m cos a — p 3 sin a /- — -r 

^4 = -<* 3 C \/ 1 + ^ Vl-sm «-/? . 

2 V mi 

Here we use the notations 5% = — sgn[sina + cos/3], 5 2 = sgn [cos (a + /?)], S 3 
sgnfcosa + sin/3] and introduce a new angle j3, 



cos /3 



m sin a + p 3 cos a 



(15) 
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It is easy to show that sin introduced in (j!4p is given by 



smp = - - . (16) 

Po - Gn 

Finally, the normalization constant of the wave function ([7]) is 



It is worth to stress that equations ((7|) , ffl~2l) , (Tl~3]) , fTH|) and ffTTl) represent 
the exact solution of the modified Dirac equation fl2]) which describes an 
electrically charged fermion (millicharged neutrino) with anomalous magnetic 
moment motion in presence of magnetized matter. 



3. Millicharged neutrino in rotating matter and magnetic field 

Consider the case of the rotating media where the matter potential / M in 
(T5]) is given by 

r = -Gn n (l,-eyu,exu,0). (18) 

Here u is an angular frequency, where e = ±1 that corresponds to parallel 
and antiparallel directions of two vectors us and B. As we shall see below, 
these two cases are quite different and may have different consequences for 
astrophysical applications. 

The Hamiltonian fll]) in the case of rotating matter becomes 

Gn 

H = 707P + 7 m + 7 S 3 /i5 2^(1 + 7s)(l + ejoliuy - e7 7 2 wx). (19) 

From the equation H^(x) = po^(x) we obtain for the wave function com- 
ponents the following equations 

(p + Gn n + p 3 )ipi - (A - 1P2 + ieGn n ure~ %ip )ip2 = (m + fiB)ip 3 , 
(p + Gn n - p 3 )tp 2 ~ (A + 1P2 - ieGn n ujre tip )i)i = (to - fiB)^, 

(Po + ^3)^3 - (A - ^A)^4 = (m + /i5)^i, 
(Po ~ Pz)^a - (A + «A)^3 = (to - fiB)tp 2 - 

It is not a trivial task to solve the system ( l20|) generally. The problem 
is reasonably simplified in the limit of small m and fiB. In fact for the 
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relativistic neutrino these parts give a negligible contribution to the energy. 
Thereby two pairs of the neutrino wave function components ipi decouple one 
from another. In this case the neutrino wave function is represented as a sum 
of two neutrino chiral states = + ^r(x), where 



1 + 75 





W 



75 



/0\ 



^3 

W 



(21) 



The second pair of equations (1201) does not contain the matter term. Thus 
right-handed chiral neutrino state, ^r, is attributed to the sterile neutrino. 
For the components "03 and ip± we obtain the solution in the following form 



V>3 
1p 4 



1 q v B 

2 V 2ttL 



P3_ M-l 



q v B 



r 2 ) e^" 1 ^, 



i quB 

2V 2ttL 



1 



q v B 



(22) 



that, in fact, describes charged particle motion in external magnetic field 
in vacuum. In this case energy spectrum is quantized due to presence of 
magnetic field and has the form 



pi + 2Nq v B, 



(23) 



where N = 0, 1, 2... correspond to the Landau levels in the magnetic field. 

The first pair of equations ( |20|) corresponds to the active left-handed 
neutrino ^l- For the components tpx and ip2 we obtain the solution in the 
form 



1 \2Gn n u - eq v B\ 
^ 1 = 2V 27Z Ji 



P3 



p^ + Gn n 



f \ZGn n uj - eqJB\ A 



Pa 



_ £ \2Gn n u - eg^E) / _ 

^ 2 ~2V 2ttL V p% + Gn n 

Here the energy spectrum is given by 



-r e 



Po = \/p1 + 2iV|2Gn n w - eg„£| - Gn n . 



(24) 
(25) 
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The solution ( I24p and ( 12 5 p describes a charged particle motion in external 
magnetic field and rotating matter. The energy spectrum (1251) represents 
the modified Landau levels that combine two different types of interaction, 
both weak and electromagnetic. Note that in the electric charge vanishing 
limit the solution ( l2"4"j) reduces to the result that was previously obtained in 



16, 26 



4. Quasi-classical treatment 

It is also interesting to consider a more general case of a millicharged neu- 
trino moving in an external environment when, in addition to the magnetic 
field, a constant electric field is presented. As it follows from (j2J) and (fl9j) . 
the neutrino energy spectrum in this case is 



p Q = ^JpI + 2N\2Gn n oo - eq v B\ + m 2 - Gn n - q<f>, (26) 

where <fr is the scalar potential of the electric field. It is possible to introduce 
the generalized "effective Lorentz force" (see also |16| ) that accounts for 
neutrino interaction with background matter and external electromagnetic 
field 

Fe// = <Ze//E e // + qeff [P X B e// ] . (27) 

Here f3 is the neutrino speed and 



Qeff^eff = \q m B m + cqB\e z 



(28) 



where g m ,B m ,E m are the matter induced "charge", "electric" and "mag- 
netic" fields correspondingly, 

q m = -G, E m = -Vn n , B m = 2n n u>, (29) 

and e z is a unit of the rotation axis. These relations directly follow from 
the exact form of the obtained energy spectrum (1261) and are in agreement 



with the most general form of the "matter induced Lorentz force" [16|. Note 
that the "effective Lorentz force" (1271) is generated by both weak and electro- 
magnetic interactions. The electric E and magnetic B fields reproduce the 
ordinary electrodynamical Lorentz force and the "matter induced Lorentz 
force" is due to E m and B m . 
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As it follows from ( 127)) , ( 125)) and (1251) , there are two types of forces acting 
on the millicharged neutrino moving in background matter with nonvanish- 
ing gradient of density in presence of external electric and magnetic fields. 
The first type is an accelerating force which is produced by the effective elec- 
tric field E e ff. The second one is the centripetal force which is produced by 
the effective magnetic field B e ^. As it follows from (127)) . the second force is 
orthogonal to the neutrino speed f3. Note that possibility of particle acceler- 



ation due to the gradient of matter density was also discussed in [16|, [29|, 130 

Consider the motion of an active left-handed neutrino that occupies one 
of the excited modified Landau levels characterized by iV ^> 1 in the case 
when the "effective electric field" E e jj ( 1251) is not presented. This situation 
may happen, for instance, for a neutrino moving in a rotating neutron star 
in the transversal to the rotating axis plane. In this case the radius of the 
neutrino quasi-classical orbit 



R= ^[r^ L dr (30) 
Jo 

depends on the parameter e and has the following form 



^ \2Gn n uj — eq u B\ ^ ^ 

It follows that for the fixed values of the magnetic field strength B and the 
rotation angular frequency u> the radius of the neutrino orbit gets different 
values for two different values of e. When the angular frequency uj is an- 
tiparallel to the magnetic field B (e = —1) the radius ([31]) decreases when 
both the angular frequency and magnetic field strength increase. Note that 
neutrinos can have bound orbits inside a rotating neutron star. The orbit 
radii of neutrino with quantum number N up to 10 10 and with energy up to 
1 eV are less than the radius of neutron star. Thus low energy neutrinos can 
be bounded inside a rotating neutron star. 

In the opposite case when two vectors oj and B are parallel (e = + 1), 
the effect of weak interactions with rotating matter and the millicharge inter- 
action with magnetic field act in opposite directions and under the condition 
2Gn n u = q u B can even "kill" each other. In this case the radius tends to 
infinity, i.e. neutrino propagates along a straight line as in a free motion. 
Note that this condition can be realized for reasonable choice of parameters 
of a neutron star (n n ~ 10 37 cm~ 3 , u = 2tt x 10 3 s _1 , and B ~ 10 12 G). 
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5. Light of millicharged neutrino (LCv) 

In more general case of rotating matter with nonconstant density and in 
presence of electric and magnetic fields a massive neutrino is moving with 
acceleration under the action of the effective Lorentz force fl28l) 



a = (GVn n + q u V<p + \2Gn n u - eq v B\ \B x eJ) — . (32) 

m 

The corresponding total power of the electromagnetic radiation due to neu- 
trino millicharge (in the quasi-classical limit) is given by 

2ql( a 2 (a/3) 2 



hcv ~ 3 V(l-/3 2 ) 2 + (l-/3 2 )V' ^ 
We tern the considered mechanism of the neutrino electromagnetic radiation 
due to neutrino millicharge, that can be emitted in the presence of back- 
ground matter and electromagnetic fields, as the "Light of (milli) Charged 
Neutrino" {LCv). 

Consider possible applications for the obtained results in different astro- 
physical environments and, in particular, in neutron stars. A neutron star 
is a very compact astrophysical object (Rns ~ 10 km) with a very dense 
neutron core and a crust around it. The crust comprises a relatively small 
fraction of total radius of the star (usually about 10%). Radial density distri- 
butions for different theoretical models of neutron stars had been considered 
in 31]. In particular, core density is extremely high (near the nuclei den- 
sity p n ~ 10 14 g/cm 3 ) and is almost uniform. In turn, crust density rapidly 
decreases from p n at the bottom to zero at the star surface. 

In case when the massive neutrino has both nonzero electric millicharge 
and magnetic moment it is interesting to compare the efficiency of the pro- 
posed new mechanizm of the electromagnetic radiation by a millicharged 



neutrino LCv with the spin light of neutrino in matter (SLv) [17|, Il8j . The 
SLv can be emitted due to the neutrino magnetic moment and was first 
proposed and studied within the quasi-classical treatment in [l9[. For neu- 
trino electromagnetic characteristics we use the values that are given by the 
present experimental limits (JTJ). 

We consider the radiation processes, LCv and SLv, during propagation 
of the neutrino from the central part of the star outwards for two specific 
regions. In the neutron core the SLv radiation power, given by 

t 4 2 2 2 77i2 1 Gn n 

I slv = -yu a m m E , a m = — , (34) 

3 ml 
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dominates over the power of the LCv. However, in the crust of the neutron 
star the rapid decrease of density is crucial. Therefore the first term of 
equation fl32|) gives a leading contribution to the LCv. From f|33|) for the 
LCv radiation power in the crust of a neutron star we obtain 



hcu = \ql{Va m ff. (35) 



Thus the ratio of radiation powers Ilcv and Islu for the neutrino mass m 
0.05 eV reads 

l LCu 



h,cu ( E 2 Vn„ 



Islu VlOd/ 3 n, 



(36) 



For the average density gradient value | CVn n | ~ 1 eV/ 1 km and the neutrino 
energy E = lOMeV we obtain 



hcu ItfeV 



Islu V Gn 



~>5 T/\ 2 



(37) 



The maximum value ~ 1 eV for the quantity Gn n can be found in the central 
core and this value decreases to the surface. 

From the estimation (|37l) it just straightforward that in case of mil- 
licharged neutrino with nonzero magnetic moment the mechanism of elec- 
tromagnetic radiation due to neutrino charge {LCv) is more effective than 
the electromagnetic radiation due to the neutrino magnetic moment (SLv) 
in dense matter, peculiar for the crust of a neutron star. 

Note that different flavor neutrinos exhibit in general different behavior in 
respect to the introduced new LCv mechanism. In fact, in the neutron star 
there is a small fraction of electrons and protons (about 10%) in addition to 
neutrons. Due to electroneutrality of matter the concentrations of electrons 
and protons are equal, n e = n p . In this case the LCv for the electron neutrino 
has the form 

hcu e = (GV(2n e - n n )f . (38) 



We expect that the direct Coulomb interaction of an extremely small neutrino 
millicharge with charged particles of matter is negligible. 

However, weak interactions of neutrino with matter with complex com- 
position depend on the neutrino flavor. And for the LCv for muon and tau 
neutrinos we obtain again (|35|) . Thus we predict different rates of LCv for 
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different neutrino ffavors. The ratio of radiation powers of neutrinos with 
the same energies is given by 

hcu e _ (Vn„ - 2Vn e )2 g 



hc^r (Vn n ) 2 

Therefore nonuniform electron component of electrically neutral matter changes 
the radiation power of the LCu e . 



6. Conclusions 

We consider electromagnetic properties of a neutrino with an electric mil- 
licharge and anomalous magnetic moment moving in an external magnetic 
field and dense matter on the basis of the exact solutions of the modified 
Dirac equation for the particle wave function in the external environment. 
Possible effects of matter rotation and effects of nonzero matter gradient are 
also included. To solve the modified Dirac equation in this background envi- 
ronment we introduce a new type of spin operator which is a superposition of 
the operators of the longitudinal and transverse polarizations. The spectra of 
the proposed spin operator and particle's energy are obtained. These exact 
solutions can be applied not only to neutrinos but to other fermions moving 
in the considered background. 

From the obtained exact solutions it follows that the neutrino energy 
spectrum is quantized due to both electromagnetic interaction of the neutrino 
millicharge with external magnetic field and weak interaction of the neutrino 
with the rotating background particles. These two phenomena overlap and, 
depending on the orientation of matter rotation axis and the direction of 
magnetic field, the effect of energy quantization can be reinforced or washed 
out. 

In the quasi-classical treatment, the energy spectrum obtained for the 
magnetized rotating matter corresponds to the neutrino motion on circular 
orbits very much similar to the Landau orbits of an electron in the external 
magnetic field. Considering background conditions peculiar for a dense ro- 
tating neutron star we predict that low energy neutrinos can be bound inside 
the star. 

We also introduce the generalized effective Lorentz force that accounts 
for the millicharged neutrino interaction with external magnetic field and 
rotating background matter with nonzero gradient of density. In general, 
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under the influence of the effective Lorentz force neutrino can be acceler- 
ated. This new possible mechanism of the electromagnetic radiation by the 
millicharged neutrino moving in magnetized nonuniform rotating matter has 
been introduced and termed the "Light of (milli) Charged Neutrino" (LCv). 
The obtained radiation power of the LCv exhibits strong dependence of the 
neutrino 7-factor (to the six order) and also on the gradient of matter density. 
It is expected that this mechanism can be efficient in different astrophysical 
environments with rather rapid decrease (increase) of matter density. 

The considered new manifestations of nontrivial neutrino electromagnetic 
properties, due to their strong dependence on the neutrino energy, are of in- 
terest for astrophysical applications, in particular in connection with the 
recently reported hints of ultra-high energy PeV neutrinos observed by Ice- 
Cube Q. 
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